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INTRODUCTION 


In this study, data is examined from two experiments 
aboard the Explorer 50 (IMP 8) spacecraft. The Johns 
Hopkins University/Applied Physics Lab Charged Particle 
Measurement Experiment (CPME) provides 10.12 second 
resolution ion and electron count rates as well as 5.5 
minute or longer averages of the same, with data sampled in 
the ecliptic plane. The Goddard Space Flight Center 
Magnetic Field Experiment provides the high time resolution 
magnetic field data (15.36 sec). IMP 8 spends about half 
its nearly circular (32 to 38 Re) orbit in the 
interplanetary field. The high time resolution of the data 
allows for an explicit, point by point, merging of the 
magnetic field and particle data and thus a close 
examination of the pre- and post-shock conditions and 
particle fluxes associated with large angle oblique shocks 
in the interplanetary field. A computer simulation has been 
developed wherein sample particle trajectories, taken from 
observed fluxes, are allowed to interact with a planar shock 
either forward or backward in time. Shock normals are 
determined from the single spacecraft method of Lepping and 
Argentiero (1) using the Imp 8 magnetic field data and OMNI 
plasma data (J.W. King, NSSDC). One event, the 1974 Day 
312 shock, is examined in detail. 

DETERMINATION OF THE SHOCK NORMAL 


The method of Lepping and Argentiero is based on the 
coplanarity theorem (2) m which the normal is determined by 
the magnetic fields: 


n = 


B,-B i )x( 


Instead of using the average values calculated from the 
data, "best estimate" values are determined xn a manner to 
be explained below. Two assumptions are necessary to use 
this method: first, that the Rankine-Hugoniot equations are 
valid for interplanetary shocks; second, that the shock is a 
step-function increase in all parameters, with superimposed 
noise. Figure 1 lends credibility to the second assumption. 
The first assumption is necessary because certain of the 
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Rankine-Hugoniot equations are used as a constraint on the 
shock normal. The equations in the shock rest frame are: 

le'Wi - « 

I 2 

(>v„v,- Sff ] =0 

(V.B. - V-.B.1* = o 

[Bi.li - « 

Eleven variables can be identified from these equations: 




Three of these variables is 
are isolated in terms of ^ J0 
the other eight. In the { 
Lepping and Argentiero *> s 
study any three variables 
could be isolated. Their 00 
data consisted of 30 sec . '* 5 
averaged magnetic field £-»<> 
data and approximately 90 i s 
sec averaged plasma data. o 
In this study 15 sec 
averaged magnetic field | o 
data was used along with 5 L 
min averaged plasma data. ->o 
Because of the large 
standard deviations 

associated with the solar 
wind data (Wx,Wy,Wz), .„ 5 
partly due to the low $ 
resolution of the data, it •* 
was necessary to isolate 
these three variables. To 10 
use them in the ^ n 

Rankine-Hugoniot equations ^ 0 
to solve for the other ** 5 
variables produced 

inconsistent results or no M 
convergence. £ , 
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The method 
eight "known" 
and the three 
calculated 


uses the 
variables 
variables 
from the 



7 9 11 1) 19 17 

H74 day 31? solsr wind 


fig 1 a) Magnetic field in solar ecliptic coordinates 
_ , __ , _ ^ „ _ __ _ Vertical lines separaata field into pra- and post- shock 

Rank me-HUgOIUOt cquacions regions which were used in the shock normal determination 
_ ^ . _ c T ^ ~ ^ b) Solar wind in solar ecliptic coordinates Shock occurs 

m a type of least squares at u 10 hours 
fit which is minimized by 

use of an iteration process. The solution should converge 
to the "best estimate" values relatively quickly. Many 
starting values should be used and these compared to assure 
convergence to the same minimum. 
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This normal determination method was tested in several 
ways. Multiple runs were made with different starting 
values (particularly, different numbers of plasma points 
were used from 15 min to 40 min of data), and these were 
checked against each other. Three events were used; one, 
the 1974 day 312 shock to be used m the remainder of this 
work, and two others which were checked against thetabn 
angles supplied by Tsurutani (private comm.) — thetabn being 
the angle between the magnetic field and the shock normal. 
The 1979 day 93 thetabn was calculated at 39° which is 
within 4° of that calculated by Tsurutani, and the 1979 day 
95 thetabn is within 6° of the Tsurutani value. The 1974 
day 312 shock was not included in the Tsurutani study. It 
was further checked by calculating the time of propagation 
of the shock from the spacecraft along the shock normal to 
the Barth and comparing against the reported SSC time. For 
a calculated thetabn of 60° and a shock speed of 470 km/sec, 
the time of arrival at Earth was 14.29 hours compared to the 
reported 14.23 hours. 

DESCRIPTION OF THE COMPUTER SIMULATION 

The computer simulation is a further development of the 
work of Chen (3). The frame of reference is as follows. A 
plane shock is assumed moving in the direction of the shock 
normal (x direction). The magnetic field is in the x-z 
plane (By=0 upstream forces By=0 downstream). The upstream 
plasma moves in the -x direction, the downstream plasma in 
the x-z plane. The electric field is in the y direction and 
constant across the shock. Data for a particular shock can 
be made to fit these conditions by suitable rotations and 
galilean transformations. The best estimate values from the 
previous section provide a consistent set for this work 
since they obey the Rankme-Hugoniot equations (as they must 
since they were determined from these equations). With this 
framework established, individual particles are allowed to 
spiral in along magnetic field lines and be reflected 
upstream or transmitted downstream due to the change m 
magnetic field at the shock surface and the initial 
conditions. The particles move according to the Lorentz 
force; each time a particle crosses the shock the final 
conditions at the crossing point become the new initial 
conditions, and the appropriate magnetic field is inserted 
into the force equation. The gain or loss in energy is 
monitored, along with the number of times the particle 
crosses the shock before being reflected or transmitted. 
Particles can be inserted either upstream or downstream from 
the shock and can move either forward or backward in time. 

DISCUSSION 

Figure 2 shows the CPME data used in the study. This is 
the 290 keV to 500 keV sectored proton data. It is highly 
structured on a short time scale particularly in _ the post 
shock region. Associated with each of the eight data 
channels is a direction consistent with the midpoint in the 



177 


SH 1.5-5 


10000 


ecliptic plane of the 45 degree sampling angle. * lasma 
collected for each of these channels will be in fc {j® 
direction. Magnetic field data is also available on this 
high resolution time scale. For several of the peaks *he 
maqnetic field is identified, and associated pitch angles 
are determined for each of the eight channels. Then, using 
a linear distribution of phase angles for each pitch, an 
ensemble of particles is sent backward in time to the shock 

surface . 

The objective of this study is a quantitative 
determination of the particle kinematics in the vicinity o£ 
interplanetary shocks, in 
particular the 1974 day 
312 shock. The 

anisotropy of the 

particle data at certain 
peak's together with the 
phase averaged energy 

changes, and comparisons, 
will be presented. 
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